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(Dated: February 7, 2008) 

The cross section for the inclusive production of isolated photons has been measured in pp collisions 
at y/s = 1.96 TeV with the D0 detector at the Fermilab Tevatron Collider. The photons span 
transverse momenta 23 to 300 GeV and have pseudorapidity \r}\ < 0.9. The cross section is compared 
with the results from two next-to-leading order perturbative QCD calculations. The theoretical 
predictions agree with the measurement within uncertainties. 
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Photons originating in the hard interaction between 
two partons are typically produced in hadron collisions 
via quark-gluon Conrpton scattering or quark-anti-quark 
annihilation 0, H H 0]- Studies of these direct pho- 
tons with large transverse momenta, p^, provide preci- 
sion tests of perturbative QCD (pQCD) as well as in- 
formation on the distribution of partons within protons, 
particularly the gluon. These data were used in global 
fits of parton distributions functions (PDFs) and com- 
plement analyses of deep inelastic scattering, Drell-Yan 
pair production, and jet production Photons from 
energetic ir and r\ mesons are the main background to 
direct photon production especially at small p7p |6| . Since 
these mesons are produced inside jets, their contribution 
can be suppressed with respect to direct photons by re- 
quiring the photon be isolated from other particles. Iso- 
lated electrons from the electroweak production of W and 
Z bosons also contribute to the background at high p'Jp. 
Previous measurements of photon production at hadron 
colliders successfully used these isolation techniques to 
extract the photon signal 0, II H El El El El 

We present, in this Letter, a measurement of the 
cross section for the inclusive production of isolated pho- 
tons with pseudorapidity |?7| < 0.9 in pp collisions at 
y/s = 1.96 TeV. (Pseudorapidity is defined as r] = 
— lntan(#/2), where 8 is the polar angle with respect 
to the proton beam direction.) The data sample corre- 
sponds to an integrated luminosity L = 326±21 nb" 1 j l4 f 
accumulated in 2002-2004 with the D0 detector |l5j at 
the Fermilab Tevatron Collider. The primary tool for 
photon detection is the central part of a liquid-argon and 
uranium calorimeter covering |?7| < 1.1. Two additional 
calorimeters, housed in separate cryostats, extend the 
coverage to < 4.2 [lif. The electromagnetic section of 
the central calorimeter (EM) is segmented longitudinally 
into four layers (EM1-EM4) of 2, 2, 7, and 10 radiation 
lengths, respectively, and transversely into cells in -q and 
azimuthal angle, A77 x A0 = 0.1 x 0.1 (0.05 x 0.05 in 
the EM3 layer at the electromagnetic shower maximum) , 
yielding a good angular resolution for photons and elec- 
trons. The calorimeter surrounds a preshower detector 
and a tracking system which consists of silicon microstrip 
and scintillating fiber trackers (0.3 radiation lengths) lo- 
cated within a 2 T solenoidal magnet. The total amount 
of material between the interaction point and the first 
active layer of the calorimeter is equivalent to approxi- 
mately 3.5 — 4.5 radiation lengths (increasing with |7y|). 
The position and width of the Z boson mass peak were 
used to determine the EM calorimeter calibration factors 
and the EM energy resolution [l7j . 

Photon candidates were formed from clusters of 
calorimeter cells within a cone of radius 1Z = 
U (A77) 2 + (A(p) 2 = 0.4; the energy was then recal- 
culated from the inner core with 1Z = 0.2. Candi- 
dates were selected if there was significant energy in 
the EM calorimeter layers (> 95%), and the probabil- 



ity to have a spatially-matched track was less than 0.1%, 
and they satisfied the isolation requirement (E to tai (0.4) — 
E EM (0.2))/E EM (0.2) < 0.10, where E total {0A) is the to- 
tal energy in a cone with TZ — 0.4 and E EM (0.2) is the 
EM energy within TZ = 0.2. Photon candidates with en- 
ergy measurements biased by calorimeter module bound- 
aries and structures were removed from consideration; 
the geometric acceptance was A — (84.2 ± 1.5)%. Po- 
tential backgrounds from cosmic rays and leptonic W 
boson decays were suppressed by requiring the missing 
transverse energy, calculated from the vector sum of the 
transverse energies of calorimeter cells, to be less than 
0.7pJ. The efficiency for the above requirements was es- 
timated with direct photons generated by PYTHIA (lif . 
Events were processed with the GEANT detector simula- 
tion package and overlaid with detector noise and min- 
imum bias interactions E3- The efficiency (excluding 
acceptance) rose from (82 ± 5)% at p\ ~ 24 GeV to 
a plateau of (92 ± 3)% at > 110 GeV. We used 
Z — > e + e~ events .17], due to the similarity between 
electron- and photon-initiated showers, to verify the se- 
lection efficiencies estimated with the Monte Carlo sim- 
ulation (MC). The photon sample was acquired with a 
three-level trigger system that relied on hardware sig- 
nals from the calorimeter and fast, software-based, pho- 
ton reconstruction. The trigger was (71 ± 9)% efficient 
for photon candidates with plj, « 24 GeV, (93 ± 2)% at 
pi rs 32 GeV and greater than 98% for p^ > 40 GeV. 
Every event was required to have a vertex, reconstructed 
with at least three tracks, within 50 cm of the nomi- 
nal center of the detector along the beam axis; the effi- 
ciency for this requirement ranged from (90.0 ± 0.3)% to 
(95.3 ± 0.1)% as a function of instantaneous luminosity. 

Four variables were used to further suppress the back- 
ground: the number of EMI cells with energy greater 
than 400 MeV within TZ < 0.2 and within 0.2 < TZ < 0.4, 
the scalar sum of the transverse momenta of tracks within 
0.05 < 1Z < 0.4, and the energy-weighted cluster width 
in the finely-segmented EM3 layer. These variables were 
input to an artificial neural network (NN) , built with the 
jetnet package E3j to suppress background and to es- 
timate the purity of the resulting photon sample. The 
NN was trained to discriminate between direct photons 
and background events. The background events, pro- 
duced with QCD and electroweak processes in PYTHIA, 
were preselected with loose criteria to increase statistics 
and to exclude high-momentum bremsstrahlung photons 
produced from partons. The resulting NN output, Onn, 
peaks at unity for signal events and at zero for back- 
ground events. Events with Onn > 0.5 were considered 
in this analysis, yielding a high photon selection efficiency 
of (93.7 ±0.2)% and good background rejection. The NN 
was tested in MC and data using electrons from Z bo- 
son decays; the resulting Onn distributions are shown 
in Fig. ^ The systematic uncertainty on the signal effi- 
ciency for the Onn requirement, estimated with electrons 
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FIG. 1: Normalized distributions of NN output (Onn) in Z 
e + e~ events for data (•) and MC (o). 



from the Z boson samples, is 2.4%. 

The photon purity (V), defined as the ratio of signal 
to signal plus background, was determined statistically 
for each p T bin. Distributions of the number of events 
as a function of Onn are shown for data and MC in 
Fig. [2 for the 44 < < 50 GeV interval. The MC 
events in this figure were weighted by the fractions that 
resulted from the fit performed with the hmcmll pack- 
age [2(J ■ The data are well described by the sum of MC 
signal and background samples, especially for events with 
Onn > 0.5. Photon purities are shown in Fig. [3] as a 
function of p T . The purity uncertainty is dominated by 
MC statistics at low p T and data statistics at high p T . 
Systematic uncertainties were estimated by using two al- 
ternate fitting functions and by varying the number of 
bins used in the HMCMLL fits. The pythia fragmenta- 
tion model was an additional source of systematic un- 
certainty. This uncertainty was estimated by varying the 
production rate of 7r°, 77, K® , and u> mesons by ±50% 
resulting in an uncertainty of 7.5% at p T « 24 GeV, 2% 
at p T « 50 GeV, and 1% for p T > 70 GeV. 

The isolated-photon cross section is measured using 
the following definition: 



NVU 



da 



dprdr] L Ap T Arj Ae 



(1) 



where N is the number of photon candidates, e is the 
combined efficiency for the selection criteria described 
above, and Ap T and Arj are the bin sizes. The factor 
U corrects the cross section for the effects of the finite 
resolution of the calorimeter. This unsmcaring was per- 
formed, as a function of p^, by iteratively fitting the 
convolution of an ansatz function with an energy reso- 
lution function. The uncertainty in this correction was 



FIG. 2: Distribution of the number of events in data (•) as 
a function of the NN output (Onn) for 44 < < 50 GeV. 
The contributions from MC background (o) and summed MC 
signal and background (□) are also shown. The MC points 
were weighted according to the fitted purity (the errors shown 
are statistical). 
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FIG. 3: Dependence of the photon purity on pip. The dashed 
line represents a fit to these points, the filled area corresponds 
to the statistical uncertainty band, and the solid lines to the 
total uncertainty band. The NN output in data was fit to the 
shapes of the MC signal and background samples. 



estimated using two different ansatz functions and in- 
cluded the uncertainty in the energy resolution. An ad- 
ditional correction was applied to p^ T for the difference 
in the energy deposited in the material upstream of the 
calorimeter between electrons (used for the energy cal- 
ibration) and photons. This correction to pip was ap- 
proximately 1.9% at 20 GeV, 1.0% at 40 GeV, and less 
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(GeV) 

FIG. 4: The inclusive cross section for the production of 
isolated photons as a function of p^,. The results from the 
NLO pQCD calculation with jetphox are shown as solid line. 



than 0.3% for pZ, > 70 GeV. The measured cross section, 
together with statistical and systematic uncertainties, is 
presented in Fig. 0] and Table [I] (The data points are 
plotted at the pt value for which a smooth function de- 
scribing the cross section is equal to the average cross 
section in the bin 22].) Sources of systematic uncer- 
tainty include luminosity (6.5%), event vertex determi- 
nation (3.6% — 5.0%), energy calibration (9.6% — 5.5%), 
the fragmentation model (7.3% — 1.0%), photon conver- 
sions (3%), and the photon purity fit uncertainty (shown 
in Fig. as well as statistical uncertainties on the de- 
termination of geometrical acceptance (1.5%), trigger ef- 
ficiency (11% - 1%), selection efficiency (5.4% - 3.8%) 
and unsmearing (1.5%). The uncertainty ranges above 
are quoted with the uncertainty at low pip first and the 
uncertainty at high p7p second. Most of these systematic 
uncertainties have large (> 80%) bin-to-bin correlations 
in p7p. Varying the choice of NN cut from 0.3 to 0.7 
changed the measured cross section by less than 5%. 

Results from a next-to- leading order (NLO) pQCD cal- 
culation (jetphox Q) are compared to our mea- 
sured cross section in Fig. 0] These results were derived 
using the CTEQ6.1M PDFs and the BFG [H frag- 
mentation functions (FFs). The renormalization, fac- 
torization, and fragmentation scales were chosen to be 
^B =^F = ^f =P T - Another NLO pQCD calculation , 
based on the small-cone approximation and utilizing dif- 
ferent FFs |28|, gave consistent results (within 4%). As 
shown in Fig.|3 the calculation agrees, within uncertain- 
ties, with the measured cross section. The scale depen- 
dence in the NLO pQCD theory, estimated by varying 
scales by factors of two, are displayed in Fig. [5] as dashed 
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FIG. 5: The ratio of the measured cross section to the theo- 
retical predictions from jetphox. The full vertical lines cor- 
respond to the overall uncertainty while the internal line indi- 
cates just the statistical uncertainty. Dashed lines represents 
the change in the cross section when varying the theoretical 
scales by factors of two. The shaded region indicates the un- 
certainty in the cross section estimated with CTEQ6.1 PDFs. 



TABLE I: The measured differential cross section for the pro- 
duction of isolated photons, averaged over \rj\ < 0.9, in bins 
of p T . (p T ) is the average p T within each bin. The columns 
<5cr s tat and Sa sys t represent the statistical and systematic un- 
certainties respectively. (Five events with p^, > 300 GeV, 
including one with pj = 442 GeV, were not considered in 
this analysis.) 



Pt 
(GeV) 


(Pi) 
(GeV) 


d 2 a 1 dp T dr/ 
(pb/GeV) 


<5o"stat 
(%) 


8<j syB t 
(%) 


23-25 


23.9 


4.14x10^ 


0.1 


23 
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26.9 


2.21xl0 2 


0.1 


19 
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31.7 


l.OlxlO 2 


0.2 


16 


34-39 


36.0 


5.37X10 1 


0.2 


15 
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41.1 
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14 
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0.4 


13 
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1.66xl0 _1 
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7.61xl0" 2 


3.5 
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3.20xl0~ 2 


5.6 
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1.59xl0~ 2 


6.5 


14 
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212 


7.36xl0~ 3 


9.8 


14 


230-300 


256 


1.81xl0" 3 


13 


15 
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lines. The span of these results is comparable to the over- 
all uncertainty in the cross section measurement. The 
filled area in Fig.|3]represents the uncertainty associated 
with the CTEQ6.1M PDFs. The central values of the 
predictions change by less than 7% when the PDFs are 
replaced by MRST2004 ^ or Alekhin2004 The 
calculation is also sensitive to the implementation of the 
isolation requirements including the hadronic fraction in 
the 1Z — 0.2 cone around the photon. The variation in the 
predicted cross section for 50% changes in the cut values 
for these criteria was found to be less than 3% [3lJ • 

In conclusion, we have measured the cross section for 
the production of isolated photons with |?7| < 0.9 pro- 
duced in pp collisions at y/s = 1.96 TeV over a wide 
range in p T , 22, < pip < 300 GeV. This extends previous 
measurements in this energy regime 0, 0, 0, 0, 
to significantly higher values of pj,. Results from NLO 
pQCD calculations agree with the measurement within 
uncertainties. 
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